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FREE EXTREME VALUES 

By Gerard Ben Arous and Dan Virgil Voiculescu 1 

New York University and University of California at Berkeley 

Free probability analogs of the basics of extreme-value theory 
are obtained, based on Ando's spectral order. This includes classifi- 
cation of freely max-stable laws and their domains of attraction, using 
"free extremal convolutions" on the distributions. These laws coincide 
with the limit laws in the classical peaks-over-threshold approach. A 
free extremal projection-valued process over a measure-space is con- 
structed, which is related to the free Poisson point process. 

1. Introduction. Free probability theory [17, 18] is a highly noncommu- 
tative parallel to a large part of basic classical probability theory. The aim of 
this paper is to add a somewhat unexpected new entry to the classical-free 
dictionary: extreme-value theory. 

Classical probability theory has constructed a complete picture for the 
behavior of extreme values of i.i.d. samples (see [7, 14]). This has been mo- 
tivated by the huge importance of this question for statistics (catastrophic 
events, insurance, reliability and, more recently, finance and risk theory). 
The results provide an exhaustive description of the possible limits of nor- 
malized extreme values as being the extreme value distribution of their do- 
main of attraction. In this note we develop a free analog to classical extreme- 
value theory. Also, surprisingly, the free extreme- value distributions turn out 
to have classical realizations already used by statisticians. 

Addition and multiplication of free noncommutative random variables 
give rise to additive and multiplicative free convolution operations on prob- 
ability measures on IR. With respect to the so-called spectral order [1] there 
are also analogs of the min and max operations for self-adjoint operators. 
Passing to the distributions, this gives rise in the case of free random vari- 
ables to the "extremal free convolution" operation that we consider here. 
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These turn out to be quite easy to express using distribution functions, 
without the analytic function machinery necessary for the other free convo- 
lutions. This can also be viewed as applying the limiting process that leads 
to idempotent analysis [8, 10] to additive free convolution. 

The natural processes in this context are noncommutative extremal pro- 
cesses over an ordered set. In particular, we show there are some remark- 
able such processes indexed by the measurable sets of a measure space: the 
projection- valued and the triangular free extremal processes over a set. We 
show that the projection-valued free extremal process over a set can be re- 
alized via the free Poisson process and therefore has a natural asymptotic 
random matrix realization. 

The analytical aspect of proving the analog of the classical results of 
Frechet, Fisher-Tippet and Gnedenko has as a main step establishing that 
the classical and free domains of attraction are actually the same. This is 
similar to the situation that occurs for the additive domains of attraction 
[3]. After finding the domains of attraction, obtaining limit laws with the 
same normalization constants as in the classical case is rather direct. The 
laws we find are called generalized Pareto laws in the statistics literature. 
More precisely, the Frechet law, the Gumbel law and the Weibull law have 
as free correspondents the Pareto law, the exponential law and the Beta law. 
Note that, reminiscent of the additive case [4], there is a certain shrinking 
of supports for the Gumbel and Weibull distributions. 

These generalized Pareto distributions, which appear as the only free 
max-stable distributions, have played a role in statistics for a long time in 
the P.O.T. (peaks-over-threshold) approach to extreme value theory [2, 12]. 

After some preliminaries on spectral order in Section 2, we discuss ex- 
tremal free convolution operations in Section 3. Then in Section 4 we in- 
troduce extremal noncommutative processes and in Section 5 we consider 
extremal free processes over a set. Section 6 is devoted to the study of the 
iteration of free extremal convolution and to the description of free max- 
stable distributions and their domain of attraction. In Section 7 we recall 
the classical P.O.T. context where the same laws appear. 

2. Preliminaries on A and V. We shall work in a tracial iy*-probability 
space (M, r), that is, M is a von Neumann algebra with an ultraweakly 
continuous faithful tracial state r. For the basic von Neumann algebra facts 
we will use, the reader may consult [9] or [16]; for the basic free probability, 
references are [17] and [18]. It will be convenient to assume M acts on 
a Hilbert space Ti. For definiteness we may take 7i = L 2 (M, r), which is 
naturally a left M-module since the L 2 -space is the completion of M viewed 
as a pre-Hilbert space with respect to the scalar product (a, b) = r(b*a). 

By Proj(M) we shall denote the set of self-adjoint projections P = P 2 = 
P* G M. If P, Q E Proj(M), the self-adjoint projections defined by (P V 
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Q)H = PH + QH and (P A Q)H =PHnQH are in Proj(M). On Proj(M) 
the order relation P < Q, that is, Q — P is a positive operator, coincides 
with PH C Q7~t, and PV Q and P AQ are the max and min with respect to 
this order. 

Let Mh = {m G M\m = m*} denote the self-adjoint elements. If a € 
and u) C R is a Borel set, then E{a\ u) G Proj(M) will denote the correspond- 
ing spectral projection. By definition, the spectral order relation [1] on 
is a -< b if E(a; [t, oo)) < £'(6; [i, oo)) for all t G R. Since a = / °° P(a; [t, oo)) 
if a > 0, it is easy to see that a -<b => a <b in case o > 0, b > and, from 
here, also for general a, b G M/j. Also a -< 6 =^ /(a) -< /(6) if / : R — > R is an 
increasing Borel function. All this also extends to self-adjoint unbounded 
operators affiliated with M. 

The operations A and V have a natural extension to (and even to 
affiliated self-adjoint operators). It is given in [1] in the case where M is 
a matrix algebra, but the extension to the case of a general von Neumann 
algebra is quite straightforward and must have occurred to many people. If 
a, be M h , then a A 6 is defined by E(aAb; [t, oo)) = E(a; [t, oo)) AE(b; [t, oo)) 
for all t G R. To see that a A b is well defined, observe that the right-hand side 
is projection-valued, decreasing and left-continuous in the strong operator 
topology as a function of t. Similarly one defines a V b by E(a V b; (t, oo)) = 
E(a;(t, oo)) V E(b; (t, oo)). Here, checking that a V 6 is well defined boils 
down to checking right continuity in t. Note that these definitions work in 
an arbitrary von Neumann algebra, that is, we did not use the tracial state r. 

In the tracial context we also have E(a V b; [t, oo)) = E(a; [t, oo)) V E(b; 
[t, oo)). Since E(c; (t — e, oo)) [ E(c; [t, oo)) as e j 0, it is immediate that the 
left-hand side is greater than or equal the right-hand side, while the fact 
that this must be an equality follows from the following inequalities that 
involving the tracial state: 

< r(E(a V b; (t - e, oo)) - E(a; [t, oo)) V E(b; [t, oo))) 

<T(E(a; (t-e,t)) + E(b;(t-e,t))). 

Using (-a) V (-6) = -(a A b) we also get E(a A b; (t, oo)) = E(a; (t, oo)) A 
E(b; (t, oo)). Passing to orthocomplements also gives in the tracial case that 

E(a A b; (-oo,t)) = E(a; (-oo,t)) V E(b; (-oo,t)), 

E(a A b; (-oo,t\) = E(a; (-oo,t\) V E(b; (-oo,t]), 

E(a V b; (-oo,i]) = P(a; (-oo,t]) A E(b; (-oo,t\), 

A(a V b; (-oo,t)) = P(a; (-oo,t)) A E(b; (-oo,t)). 

Note, that since only spectral projections are involved, the definitions and 
properties of A and V also extend to unbounded self-adjoint operators affili- 
ated with (M, r). Note also that if / : R — > R is an increasing Borel function, 
then /(a) A f(b) = f(a A b) and f(a) A f(b) = f(a V b). 
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Lemma 2.1. If P,Q G Proj(M) are freely independent in (M, r), t/ien 
r(PvQ) =min(r(P) + r(Q),l) and r(PAQ) = max(0, t(P) + t(Q) - 1). 

This is not a new result. It can be obtained using free convolution, addi- 
tive or multiplicative. Indeed, P A Q = E(P + Q, {2}) = E(PQP, {1}). Since 
the Cauchy transforms of the distributions of PQP and P + Q are alge- 
braic, r(P A Q) is given by the residue at 1 or, respectively, at 2 of the 
corresponding Cauchy transform (see 3.4.1 in [17] or 3.6.7 in [18]). 

The free independence of two projections of given trace completely deter- 
mines the trace on the algebra they generate. Since in this paper we focus on 
P A Q and P V Q, it will be convenient to also have at hand a more relaxed 
concept. 

Definition 2.2. Two projections P,Q G Proj(M) are in general posi- 
tion if the equivalent conditions r(P V Q) = min(r(P) + r(Q), 1) and r(P A 
Q) = max(0, r(P) + r(Q) — 1) are satisfied. We will also say that two un- 
bounded self-adjoint operators a, b affiliated with M are in general spectral 
position if, for each t G ffi, the projections E(a; [t, oo)) and E(b; [t,oo)) are 
in general position. 

The preceding definition for self-adjoint operators may seem to depend 
on choosing [t, oo) instead of (i,oo), but this is actually inessential. 

LEMMA 2.3. If the unbounded self-adjoint operators a,b affiliated with 
(M,t) are in general position, then E(a; (t,oo)) and E(b; (t, oo)) are in gen- 
eral position. In particular, —a and —b are in general position. 

The preceding lemma follows from the inequalities 
< T(E(a; (t, oo)) V E(b; (t, oo)) - E(a; [t + e, oo)) V E(b; [t + s, oo))) 
< r(P(a; (t, t + e)) + E(b; (t, t + e))). 
Returning to the spectral order we have the following lemma. 

Lemma 2.4. Ifa,beM, a > 0, b > 0, then (2~ l {aP + &p)) 1 /p | «V b as 
p — ► +oo. 

With the superfluous condition that M be finite-dimensional, this is Lem- 
ma 6.15 in [1]. The proof in [1] that 2 -1 / p (a p + b p ) l / p is increasing, works 
in general and the rest of the proof also works with minor adjustments. 
If X denotes the limit, then 2~ 1 (a p + bf) < (a V b) p gives X < a V b since 
t — s- is operator- increasing for p > 1. On the other hand, if k > 1, X k = 
s - limp^ 2- k / p (a p + 6P) fc / p = s - limp^ 2~ 1 / p (a fcp + b kp ) x l p > a k , so that 
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if ££E(X; (-00, t])H, then t k \\£\\ 2 > (X k £,£) > (a k £,£) for all k G N, so that 
£ € E (a; (- 00, t])H. Similarly, £ G E(b; (-00, t])H, so that E(X;(-oo,t]) < 
E(a; (—00, t]) A E(b; (—00, t\). This easily gives X > aV b. 

Corollary 2.5. If a,b e Mh, then 

a V b = s — lim p _1 log(exp(pa) + exp(p6)). 

— >oo 

Indeed this follows from log(exp(a V b)) = a V b, 2~ 1 ^ p —* 1 and the pre- 
ceding lemma applied to expa, expfe. 

3. Extremal free convolutions. By Prob(R) we shall denote the proba- 
bility measures on R and by Prob c (R) denote those with compact support. 
If /t G Prob(R), then F(t) = //((— oo,£]) is its distribution function. 

Definition 3.1. If G Prob(R) have distribution functions F, G, then 
the upper and the lower extremal free convolutions y,Mv and, respectively, 
//Eli/ are given by the distribution functions H(t) = max(0,F(t) + G(t) — 1) 
and, respectively, K(t) = min(F(i) + G(t), 1). 

Proposition 3.2. Ifa,b are self-adjoint unbounded operators affiliated 
with (M, r) that are in general spectral position and if /t a , /*& G Prob(R) are 
their distributions, then /t a M/i6 and /t a EI/i& are the distributions of aVb 
and, respectively, a A b. 

The preceding proposition is immediate from the definitions. Using Lem- 
ma 2.1, Definition 2.2 and the adaptation of free independence to unbounded 
operators affiliated with (M, r) [5] we have the following consequence. 

Corollary 3.3. If a,b are freely independent unbounded self-adjoint 
operators affiliated with (M,t), then the distributions ofaVb and aAb are 
equal to // El/if, and, respectively, to fj, a \Sfib where /i a) A*b are the distributions 
of a, b. 

The extremal free convolution operations also have natural descriptions 
without invoking the distribution functions. 

Proposition 3.4. Let fi,u G Prob(M) and let t = inf{x G R|(/t + v) x 
((x, 00)) < 1} and s = sup{y G R|(/i + i/)((— 00, y)) < 1}. Then 

fj,Eftu = X(t,oo) ■ i^ + v) + (1 - + z/)((i,oo)))<5 4 

and 

/i0z/ = X(_ OO)S ) • + i/) + (1 - + z/)((-oo,t)))5 s . 
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Proof. Given that the map x — > —x interchanges the two operations, 
it will suffice to check the first equality. If F, G and H are the distribution 
functions of fi, v and the right-hand side of the equality we are checking, 
respectively, it is immediate that 

1 - H(x) = min((l - F{x)) + (1 - G(x)), 1), 

which gives 

H(x) = max(F(x) + G(x) -1,0). □ 

Translating Corollary 2.5 into free convolutions language, we obtain the 
following result. 

Proposition 3.5. Let /j,, u e Prob c (M). Then fiEAis is the weak limit of 
(fc _1 log)*((exp£;«)*/x EB (expfc»)*i/) 

as k — > +oo. 

Note that the preceding proposition can be given additional precision 
using Lemma 2.4 with the 2 _1//p factor. 

Remark 3.6. It is interesting to note that in the classical context [7, 
14], the operation on probability measures, analogous to 0, corresponds to 
the multiplication of the distribution functions. Thus, viewing distribution 
functions as functions with values in [0, 1], the passage from classical to free 
probability means replacing the multiplicative semigroup ([0,1], •) by the 
semigroup on [0, 1] that results from the binary operation (s,t) — > max(0, s + 
t-1). 

4. Extremal noncommutative processes indexed by an ordered set. 

Definition 4.1. Let (J, <) be an ordered set and let (M, r) be a tracial 
iy*-probability space. An M-valued upper extremal process indexed by J 
is a family (Y(i))i^j of self-adjoint unbounded operators affiliated with M 
such that if a 6 J is the least upper bound of a set uj C J ', then 

Y{a) = V{Y(i)\i £w) 

[for sets, V is defined like for pairs on the spectral projections E{-\ (t, oo))]. 
The process is projection- valued if Y(i) £ Proj(M), iEj. 

Remark 4.2. The extension of the supremum V to infinite families no 
longer has the property that spectral projections for (t, co) may be replaced 
by spectral projections for [i,oo). For instance, \/ neN (l — l/n)I = I, but 
E({1 - l/n)I; [1, oo)) = while E(I; [1, oo)) = /. 
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Remark 4.3. The preceding definition implies that if a,j3 G S, ct<j3, 
then Y(a) -< Y{0) with respect to the spectral order. Indeed, it suffices to 
apply the definition to to = {a,/3}. 

If Y\, Y2, are M- valued upper extremal processes indexed by the same 
ordered set J', then we may define Y\ V Y2 by (Yi V Y2)(a) = Y±(a) V ^(a), 
which is also an upper extremal process indexed by J . Clearly the construc- 
tion extends to families of processes. 

If Yj, j = 1,2, are upper extremal processes indexed by J with values in 
von Neumann algebras (Mj,Tj), j = 1,2, we may construct their free sup 
as the process ci(Yi) V 172(^2)) where aj :Mj — > Mi * M2 are the canonical 
inclusions into (Mi * M2,ti * T2). Again, more generally we may take the 
free sup of a family of processes. 

Consistent with these considerations is the spectral order on processes 
Yi -< Y 2 defined as Yi(a) -< Y 2 (a) for all aej. 

Remark 4.4. If Y is an M-valued upper extremal process indexed by 
then j9a-> E(Y(a); (t,oo)) G Proj(M) is a projection-valued upper 
extremal process indexed by J . 

Example 4.5 (Spectral measures). If a is an unbounded self-adjoint 
operator affiliated with (M, r), then 

R 3 t -> E(a; (-00, £)) G Proj(M) 

is a projection-valued upper extremal process indexed by M. 

Remark 4.6. Note that a -< b is equivalent to E(a; (— 00, t)) < E(b; 
(— 00, t)), t G E. Thus, Ando's definition of the spectral order [1] can be 
interpreted as replacing the operator by the projection- valued process of 
Example 4.5 and using the natural order on such processes derived from the 
order on projections. 

Remark 4.7. The definition of upper extremal processes can be ex- 
tended to processes where the Y(a) are self-adjoint unbounded operators 
affiliated with M. Note also that if / : M. — s> R is an increasing function that 
is lower semicontinuous, then 

j3a^f(Y(a)) 

transforms an upper extremal process into another upper extremal process. 
We denote this process by f*Y. 
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5. Free extremal processes over a set. Let (X,B,fM) be a measure space 
with \i a positive cr-finite measure. The free 'projection-valued upper extremal 
process over (X,B,/i) will be a projection-valued upper extremal process 
B/^3 lo — > Y(lo), where £>/~ denotes B modulo null sets with values in 
some (M, r) so that if {uJk}keK are pairwise disjoint (modulo null sets), 
then the y(wfc), k € K, are freely independent and t(Y(lu)) = min(^(u;), 1). 

It is easy to see that it is sufficient to construct such a process when 
X is a finite set; the general case then is roughly the result of viewing 
X as arising from an inverse limit using finite partitions, while M would 
be constructed as a direct limit (we leave the details to the reader). If X 
is finite and X = {x\, . . . , x n }, {xj} € B, then let P±, . . . , P n & Proj(M) be 
freely independent with r(Pj) = min(/i({xj}), 1). We then define, for uj C X, 

{j\xj€u} 

That t(Y(uj)) = mm(fi(uj), 1) is then a consequence of Lemma 2.1. 

It is easily seen that under the additional requirement that M be gen- 
erated by {y(u;)|u; € B}, the free projection-valued upper extremal process 
over (X,B,fJ,) is unique up to isomorphism. 

The free projection-valued upper extremal process over (X,B, //) is re- 
lated to the free Poisson process over (X,B,fj,). To explain this, we start by 
recalling some facts about the free Poisson process. Based on realizations 
of free Poisson variables using semicircular or circular elements (see [11], 
Remark 1.7, Corollary 1.8 in the main text and Lemma 1.4, Remark 1.5, 
Theorem 1.6 in the Appendix), there are constructions of free Poisson pro- 
cesses over a set (see Section 6.2 and 1° in Section 6.6 of [17]). We will use 
the construction that involves circular elements, which we summarize in the 
following theorem. 

Theorem 5.1 ([11, 17]). Let (C t ) te / be *-freely independent circular 
variables and let (OJ tg / be spaces of events with a-algebras S t and proba- 
bility measures v L . Let further £,//) be the disjoint union of E t , v^). 
Let (A l ,t l ) denote the W* -probability space L°°(£l L , E t , u L ) with r t the expec- 
tation functional given by v L . Assume (A L ) ie j and ({C t , C*}) L( zj are freely 
independent and contained in (M, r). If a 6 S, v{a) < oo, let 

n(a)=£c tXa na.c;, 

where Xanfi t G A L is the indicator function of a n Then: 

(i) 11(a) has free Poisson distribution with parameters a = v{a), 6=1 
in (M,r). 

(ii) If a is the disjoint union of € X, v{a) < oo, then 11(a) = SfceNH(afc)- 
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(iii) If ak G X, v{otk) < oo are pairwise disjoint, k G N, then (U(ak))k<=N 
are freely independent. 

Using the facts on freeness of Gaussian and deterministic diagonal ma- 
trices, it was also noted (see [17], the end of Section 7.3) that one obtains 
asymptotic random matrix realizations of the free Poisson processes. 

To avoid complicating notation, we gave in Theorem 5.1 a construction of 
the process for f2, which can be represented as a disjoint union of probability 
measure spaces. Clearly if (X, B, fi) is a u-finite measure space, we can realize 
the free Poisson process over {X, B, fi) by finding f2 as above so that f2 D X, 
BcH and v\X extends \x. 

Let (X,B,fi) be a cr-finite measure space and let II: Bf — ► M^, where 
Bf are the sets in B with finite measure, be the free Poisson process, that 
is, 11(a) is a free Poisson noncommutative random variable with parame- 
ters n(a) and 1, and for disjoint a±,a2, ... so that Ylkflip-k) < oo, we have 
J2k^( a k) = n(Ufc a fe) an d the n(o!fc) are freely independent in (M, r). The 
following theorem gives the connection between the two kinds of processes. 

Theorem 5.2. IfY{a) denotes the range projection ofH(a), then Y{a) 
is the free projection-valued upper extremal process over (Af, £>,//). More- 
over, if ji has no atoms, then the von Neumann algebras {H{a)\a G B} and 
{Y(a)\a G B} are equal. 

Proof. In view of the formula for the distribution of a free Poisson ran- 
dom variable (Section 2.7 in [17]), we have t(Y(o)) = min(/i(a), 1). Also, 
clearly 11(a) depends on a only up to null sets and if ai,a2,... are dis- 
joint, then Y{a\),Y{a.%), . . . are freely independent, since H(ai),H(a2), ■ ■ ■ 
are freely independent. Further, if a\, ct2, ■ ■ ■ are disjoint and a = [j k a^ 6 Bf, 
then Y(a>) = \J k Y(ak)- Indeed, Y(a) is the projection onto (kern(a))- 1 - and 
it suffices to show that ker 11(a) = /\ fc kerLI(afc), that is, = nfckerll(afc). We 
have 

£ G ker LT(a) £ G ker Il(a) 1/2 

O (II(a)£,£} = 0^]T(Il(a fc )£,£} = 

k 

^(II(a fc )£,0=0, fcGN 

& £ G flkern(afe) 1 / 2 = f|kerll(a fc ). 

k k 

Our use of N as an index set is no loss of generality since 11(a) depends 
on a only up to null sets and /U is <7-finite. Thus we have checked that Y is 
the free projection-valued upper extremal process. 
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In view of the formula for the free Poisson distribution (Section 2.7 in 
[17]) we have that if /x(a) = a < 1, the spectrum of 11(a) is {0} U [(1 — 
a 1 ! 2 ) 2 , {I + a 1 / 2 ) 2 } and Y(a) is the spectral projection of 11(a) for [(1 — 
a l l 2 ) 2 ,{l + a 1 ' 2 ) 2 ]. This gives \\Y(a) - U(a)\\ < (1 + a 1 / 2 ) 2 - 1 < 3a 1 / 2 and 
\Y(a) - n(a)|i < \\Y(a) - U(a)\\T(Y(a)) < 3a 1 / 2 • a = 3a 3 / 2 . 

To prove that the von Neumann algebras generated by {y(a)|a G B} and 
{n(a)ja G £>} coincide, it suffices to show their L 1 -spaces coincide, and hence 
it suffices to show 11(a) is in the L 1 -closure of the linear span of the Y(a). 
Given a G Bf, if fi is diffuse, there are pairwise disjoint a±, . . . ,ak G Bf so 
that ai U • • • U afc = a and /J>(ctj) < e. Then we have 

[1%) - (Y( ai ) + ■■■ + F(a fc ))li < E l n ( a i) " 

3 

<3e 1 / 2 ( / u(a 1 ) + --- + / u(a fc )) 
< 3e 1/2 ^(a), 

which proves our assertion. □ 

Remark 5.3. Note that in the realization of the free Poisson process 
via circular elements, Y{ot) being the range projection of Sc t x(a H £\)c* 
is the same as being the range projection of Cx(a), where C is a column 
matrix with entries C L and x( a ) is a diagonal matrix with entries xi a ^ ^t) 
on the diagonal. This can also be translated into the asymptotic random 
matrix realization (see Section 7.3 in [17]): the Y(a) are the large iV limit 
of rArDjv(a) [or equivalently of rAri?Ar(a)r^-], where T^r is an appropriate 
N x [tN] matrix with i.i.d. Gaussian entries and D^{a) are deterministic 
[Nt] x [Nt] projection matrices with joint limit distribution the same as the 
x(a)'s in L°°(f2, S, t~ l v) [here t = < oo; we leave to the reader the easy 
adaptation to = oo]. 

Remark 5.4. By relaxing the condition Y(u) G Proj(M) to Y(u) G 
Mh, we may consider more general free upper extremal process over a 
set, that is, in addition to the condition that uj — > Y(u>) be an upper ex- 
tremal process over we also require that, for pairwise disjoint {uJk}keK, 
the (Y(ujk))keK be freely independent. Note that in view of Remark 4.4, 
E(Y(u>); (t, oo)) will then be a projection- valued upper extremal process in- 
dexed by and the freeness requirement for disjoint {uJk\k&K also carries 
over to the E(Y(u>k); (t, oo)). 

An example of such a more general free upper extremal process over 
(X,B,fj,) is the free triangular upper extremal process Z(cj) over (X,B, fj,). It 
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is characterized by the fact that E(Z(u); (t, oo)) = Rtioj) has the properties 
t <0^R t (u) =1, 
t>l=>Rt(u})=0, 
< t < 1 => T(R t (u)) = min((l - t)n(u), 1). 

The kind of argument that we used for the existence of the projection- 
valued free process over (X,B,fj>) also works for the triangular process (we 
leave the easy details to the reader). 

Note also that, in addition to the triangular process, free upper extremal 
processes over (X,B,fj.) can be obtained applying Remark 4.7 to Z(u), that 
is, processes f*Z, where / is an increasing lower semicontinuous function 
R^R. 



6. Free max-stable distributions and free max-domains of attraction. It 

will be convenient to work with distribution functions and to adapt some of 
the notation and definitions of Section 3. 



Definition 6.1. If F and G are two distribution functions on the real 
line, we define the distribution function FEUG to be (F + G — 1)+ and we 
define the n-fold iterate of this operation to be 

F Un = FM_-MF = (nF - (n - 1)) + . 

n times 

We want to study the possible asymptotic behavior of F an when n tends 
to infinity. Let us start with some trivial properties of the tail and the 
support of F Sn . If F = 1 — F denotes the tail of the distribution function 
F, then 



FMG=(F + G) Al 

and 



F Un = nF A 1. 

If [a(F),u)(F)] denotes the support in [—00,00] of the probability distribu- 
tion defined by F, then for all n > 2 and F, 

a(F an ) > -00 

and 

oo(F EZn )=Lo(F). 

More precisely, 

a(F an ) = sup jx G R\F(x) < 1 - - j 
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and 

lim a(F mn ) =u(F). 

n—>oo 

In case there is u n such that 

F(u n ) = l-n-\ 

then 

a{F® n )>u n 

and one can easily interpret F an as a conditioned (or thresholded) distri- 
bution. Indeed, if F is the distribution function of a random variable X , 
then 

F Un {x) =P(X <x\X>u n ), 

that is, F® n is the distribution function of the random variable X condi- 
tioned to be larger than the threshold u n (cf. Proposition 3.4). 

Definition 6.2. A distribution function F is freely max-stable iff for 
every n > 1, there exist a n ,b n € M, a n > 0, such that 

F Mn (a n x + b n ) = F{x). 

Remark 6.3. If F is freely max-stable, then its support is bounded 
from below, since a(i ?02 ) > — oo. 

Definition 6.4. A distribution function F is in the free max-domain 
of attraction of the distribution function G if there exist a n ,b n € M., a n > 0, 
such that, as n — > oo, 

F > n x + & n )^G(*) 

(i.e., convergence at every point of continuity of G). The free max-domain 
of attraction of G will be denoted by Domf ree (G'). 

Theorem 6.5. The following statements are equivalent: 

(i) G is freely max-stable; 

(ii) Dom frcc (G) ^ 0; 

(iii) G G Domf rcc (G). 

Proof. Clearly (i) => (iii) (ii) and we are left with proving that (ii) 
(i). If F G Domf ree (G), then there are a n , b n € R, a n > 0, so that at every 
continuity point x of G, we have 

lim nF{a n x + b n ) A 1 = G(x). 
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Hence, if k > 1 is an integer, then 

lim nF(a n x + b n ) A AT 1 = G(x) A AT 1 

n — >oo 

and 

lim nF(a nk x + b nk ) A A;" 1 = k~ 1 G(x). 

Thus, by a slight generalization of Khintchine's law of types, which we give 
in Lemma 6.6, we infer the existence of a k ,(3 k £R, ct k > 0, such that 

G(a k x + h ) A AT 1 = k~ l G{x) 

or, equivalently, 

kG(a k x + j3 k ) Al = G(x), 

which is the same as 

G m (a k x + (3 k ) = G(x). 

So, modulo Lemma 6.6 we have proved G is max-stable. 

To formulate the extension of Khintchine's law of types, let us call G a 
c-defective distribution function (here < c < 1) if G is a nondecreasing, 
right continuous function on R so that 

lim G(x) = 1 

x— >+oo 

and 

lim G(x) = c. 

X— >— oo 

Then G can be viewed a the distribution function of a random variable that 
takes the value — oo with probability c. 

Lemma 6.6. If F n is a sequence of c-defective distribution functions and 
if G and G> are c-defective, nondegenerate distribution functions such that, 
as n — ► oo, we have 

F n (a n x + b n ) ^>G(x), 
F n (a n x + (3 n )-^ G*(x) 
for real constants a n > 0, b n ,a n > 0, (3 n , then 

lim — = a and lim = b 

for some a > 0, b € R and 

=G(ax + 6). 
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This lemma is a trivial consequence of Khintchine's law of types. Indeed 
one reduces it to the usual result by defining 

~ G{x)-c ~ G*(x)-c 

G(x) = — and GJx) = — . 

1 — c 1 — c 

To conclude the proof of the theorem, we apply the lemma to the se- 
quence of (1 — /c -1 )-defective distribution functions 1 — (nF A k~ l ) and to 
the (1 — /c _1 )-defective distribution functions 1 — G A k~ l and 1 — Gk . □ 

Definition 6.7. We will say that a distribution function F is of free 
extreme-value type if F is of the same type as one of the following classes 
of distributions: 

Type I: The exponential distribution F(x) = (1 — e -x )+. 

Type II: The Pareto distribution F(x) = (1 — for some a > 0. 

Type III: The Beta law F(x) = 1 — \x\ a for —1 < x < and some a > 0. 

Theorem 6.8. G is freely max-stable iff G is of free extreme-value type 
[i.e., there exist a > and b real constants so that G(ax + b) is one of the 
distributions listed in Definition 6.7]. 

Proof. Each of the distributions in Definition 6.7 is freely max-stable 
by a straightforward computation. To prove the converse statement, we need 
the following lemma. 

Lemma 6.9. If G is freely max-stable, then there exist measurable func- 
tions a(s) > 0, b(s), where s £ [1, oo), so that 

G as {a(s)x + b(s)) = G(x) 

with G as defined for s 6 [1, oo) by 

G® s (y) = sG(y) Al. 

Proof. Whereas G is freely max-stable, there exist real numbers a n > 
and b n such that 

G(x) = nG(a n x + b n ) A 1. 
If u s (x) = a [ns] x + b [ns] , then 

G{x) = [ns]G{u s {x)) Al. 

This easily gives 

lim nG(u s (x)) A s" 1 = G(x)s~ 1 . 
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Thus, we have 

nG(a n x + 6 re ) A s" 1 = G(x) A s" 1 

and 

lim nG(a[ ns p + b[ ns j) A s" 1 = G'(x)s _1 . 

By the extension of Khintchine's law of types in Lemma 6.6, we conclude 
that there exist measurable functions a(s) > and b(s) such that 

s~ 1 G{a(s)x + 6(a)) = G{x) A s' 1 

or, equivalently, 

G(a(s)x + b(s)) = sG(x) A 1, 

that is, 

G(o(s)x + 6(s)) = G 0s (2;). □ 

By applying Lemma 6.9 to G® st , one easily gets that 

a(st) = a(t)a(s), 
b(st) = a(t)b(s) + b(t) 

if s,t € (1, oo). 

It is easily seen that if one extends a(s) to (0, oo) by a(l) = 1 and a(s) = 
(a(s -1 )) -1 if s < 1, then the equation a(st) = a(t)a(s) holds for all s > and 
t > 0. Since a(s) is measurable, it is known that there exists ^GR such that 
a(s) = s 9 . 

Case I. If 9 = 0, that is, a(s) = 1, then 6(st) = 6(s) + b(t) for s,i > 1 and 
6(s) = e 6 ^ s ) satisfies 6(si) = b(s)b(t). Thus again, there exists c € M such that 
b(s) = — clns and 

G 0s (x) =G(x-lns) 

for s > 1. It is easy to see that this implies that G is of the same type as the 
exponential distribution. Indeed, we just check that c > 0, since 

G{x - cln s) = W s {x) = sG(x) A 1 > G(x), 

which implies that x — clns < x and thus c > 0. If G is nondegenerate, we 
cannot have c = 0. Then we see that G(x) < 1 for all x. Indeed, if G(x) = 1, 
then G{x) = 0, so that GF°(x) = 0; that is, G Us {x) = 1 for all s > 1. Hence 
G(x — clns) = 1 for all s > 1, which gives G(y) = 1 for all y < x, which is 
not possible. 

Now, if a = G(x) > and y = x — clns, then 

G(y) = sa A 1 = (ae x/c e- y/c ) A 1 
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so that a{G) = x + clna. This shows that 

lnG(x) = (a(G) -x)c _1 for x > a{G), 

that is 

G(z) = e ^{G)-x)c-^ 

and G(a(G) + ex) is exponentially distributed. 

Case II. If 9 > 0, then a(s) = s e > 1 for s > 1. From 

6(si) = a(t)b(s) + 6(t) = a(s)b(t) + &(s) 

we see that b(s)(l — a(s))" 1 is a constant c. We infer that b(s) = c(l — s _e ) 
and G 0s (x) = G(s- fl (a; - c) + c). If #(x) = G(x + c), then 

# 0s (x) = #(s-V), 

that is, sH(x) A 1 = H{s~ e x). If x < is such that < H(x ) < 1, then 
this shows that H(s~ e xo) is an increasing function on some interval (1, 1 + 
e), since it is equal to sH(xq) for s close to 1. However, this function is 
nonincreasing if xq < 0. Thus for every x$ < we see that H(xq) € {0, 1}. 
Since H is nondegenerate, it is then impossible that H (0) = 0. Choosing 
s > 1 such that s _1 < H(0), we have 

1 = sH(0) A 1 = H(s~ e • 0) = tf(0). 

Hence ^(0) = 1. This proves that a(H) > 0. 

Let xo > be such that < H{x$) < 1. Then if y = s~ e xo, one sees 
that H(y) = sH(x ) A 1. Hence = 1 iff s > s = (H(xo))- 1 , that is, 

iff y < So ^o- This proves that a(H) > and a(H) = Sq®xo. This is valid for 
every xo such that < H(xq) < 1, so that a(H) = (H(x)) e x for x > a(H) 
and 

H(x) = (aiH)) 1 /^- 1 / 6 , 

that is, H is a Pareto distribution. 

Case III. If 9 < 0, this is completely similar to Case II and G is of the 
type of a Beta law. □ 

In particular we have also proved the following fact: 

Theorem 6.10 (Free extremal type theorem). The following statements 
are equivalent: 

(i) There exists a distribution F and constants a n ,b n € K, a n > 0, such 
that F Mn {a n x + b n ) ^ G{x) asn^oo. 

(ii) Distribution G is of the type of a free extreme-value distribution. 
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This result in the free probability setting is the equivalent of the classical 
extremal type theorem. One can introduce a natural mapping that relates 
the two statements. 

Given c > 0, we define a function on [0, 1] by 

f c (u) = (l + clnn)+. 

Then f c is nondecreasing and f c (u) = iff u < e _1 / c . Moreover / c (l) = 1. 
We have 

f c (uv) = (1 + chili + clnt;) + 

= (fc(u)+f c (v)-l) + . 

We can endow the set [0, 1] with two semigroup structures: one that arises 
from usual multiplication and the other that arises from the operation re- 
lated to free convolution: 

u*v = (u + v — 1) + . 

Then f c (for any c > 0) is a homomorphism between these two semigroups 
on (0, 1] . 

This homomorphism gives rise to a homomorphism between the semi- 
groups of probability distribution functions endowed with either pointwise 
multiplication or the operation M, which coincides with performing the op- 
eration * pointwise. It follows that 

f c (FG) = f c (F)Mf c (G) 

so that 

f c (F n ) = f c (Ff n . 

It is clear that the free extreme-value distribution functions are obtained 
from the classical ones by the map f\. It is also clear that if F is classically 
max-stable, then f c (F) is freely max-stable. If one could prove the converse 
directly, that is, that if G is freely max-stable, it is the image by f c of a 
max-stable distribution-function F, then we could derive the free extremal 
type theorem directly from the classical one. We were able to obtain this 
result via a different route. 

We will now prove that free max-domains of attraction and classical max- 
domains of attraction coincide for corresponding laws and that the normal- 
izing constants are also equal. 

We begin with Type I. 

Theorem 6.11. The free max-domain of attraction of the exponential 
distribution coincides with the classical max-domain of attraction of the 
Gumbel distribution and the normalizing constants are equal. 
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Proof. F is in the free max-domain of attraction of the exponential 
distribution iff there exist b{s) and a(s) > for s > 1 such that for all x > 0, 
we have 

lim sF(b(s) + xa(s)) = e~ x . 

s—>+oc 

(This is clear for s running over N and we can take a(s) = a([s]), b(s) = 
b([s]).) In particular, lim s _> +0O sF(b(s)) = 1. Hence we may assume b(s) < 
uj(F) and that b(s) is nondecreasing. It is easy to see that this implies that F 
is in the domain of attraction of the exponential distribution iff there exists 
a function g(t) > such that, for all x > 0, we have 

lim nt + * 9 {t)) =e _ x 

tTw(F) F(t) 

If t < uj(F), let us denote by U the function U(t) = l/F(t) so that 

lim ^ + ft)) =e x 
f\w{F) U(t) 

and U is T-varying (see [14], the definition in 0.4.3, page 26) on (a(F),uj(F)). 
By [14], Proposition 0.10, page 28, this is equivalent to F being in the 
classical max-domain of attraction of the Gumbel law A. □ 

We now turn to the free max-domains of attraction of the Pareto distri- 
bution. 

Theorem 6.12. The following statements are equivalent for a > 0: 

(i) F is in the free max-domain of attraction of the Pareto distribution 
with exponent a. 

(ii) F is in the max-domain of attraction of the Frechet distribution $ a . 

(iii) F is —a-regularly varying at oo. 

Moreover the normalization constants in (i) and (ii) can be chosen of the 
form a n = u n , b n = 0. 

Proof. The equivalence of (ii) and (iii) is a classical fact due to Gne- 
denko (see [14], the definition in 0.4.1 on page 13 and Proposition 1.11 on 
page 54). 

To prove (iii) =>■ (i), assume F is —a-regularly varying at oo, that is, 

lim -^—- L = x 

t->oc F(t) 

for all x > 0. Then u n = ml{t € R\F(t) < n" 1 } will be such that 

lim nF(u n ) = 1 
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and hence 

lim nF(u n x) = x~ a 

n~ >oo 

for all x > 0. It follows that 

lim F Un (u n x) = lim F(u n x) A 1 = x~ Q A 1 

n—>oo n— >oo 

if x > 0, and since F® n (u n x) is decreasing and less than or equal to 1, the 
limit will be 1 if x 6 (— oo, 1). This proves (i). 

To conclude the proof we will show that (i) =4> (iii) . If for some choice of 
constants a n , b n € R, a n > 0, we have 

lim F^"(a n x + 6 n ) = x~ a A 1 

n^oo 

for x > 0, then if J7= 1/F, 

lim n~ x U(a n x + = x a when cc > 1. 

n— *oo 

Then if V(y) = U*~(y) = inf{s :U(s) > y}, we have 

lim^-fc^a-^y 1 / ify>l. 

n — >qo 

With a(t) = aw and 5(£) = &w, we then have 

lim(F(ty)-^))a(t)- 1 =y 1 /« 

I— >oo 

and 

lim - y(ty 2 ))a(t)" 1 = - y l l a 

t— >oo 

when y,y 1} y 2 > 1. 
This implies 

lim a(ix)a(£) -1 

t— >oo 

= lim ((F(txyi) - y(tey 2 ))a(t)- 1 )/(y(tey 1 ) - F(txy 2 )a(te)- 1 ) 

I— *oo 

= ((xyx) 1 /" - ( X y 2 )^)/(y{ /a - y\ la ) = x 1 ^. 

Then, for any fixed y\ > y 2 > 1, the function V(ty±) — V(ty2) is 1/a- 
regularly varying at oo as a function of t. Since a function Wit) is 1/a- 
regularly varying at oo iff W{tz) for some fixed z > is 1/a-regularly varying 
at oo as a function of t, we infer that V(ty) — V(t) is 1/a-regularly varying 
at oo as a function of t for all y > 0, y 7^ 1. Thus we can use the last part 
of the proof of Proposition 1.11 in 1.2 of [14] (starting with the last two 
paragraphs of page 55 and continuing on pages 56 and 57) and conclude 
that F is regularly varying of index —a at 00. □ 

We finally turn to the domain of attraction of the Beta law (Type III). 
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Theorem 6.13. The following statements are equivalent: 

(i) F is in the free max-domain of attraction of G(x) = 1 — \x\ a for 
-1 < x < and a > 0. 

(ii) F is in the classical max-domain of attraction of the Weibull distri- 
bution fy a . 

(iii) uj(F) < oo and F(u(F) — u) is regularly varying of exponent a at 0, 
that is, lim^jo F p(u(F)-h) = xa if x > 0. Moreover, the normalization con- 
stants can be chosen to be 

a n = uj(F)—u n and b n = u>(F). 

Proof. The equivalence of (ii) and (iii) is a classical fact due to Gne- 
denko (see [14], Proposition 1.13 in 1.3). 

To prove that (iii) implies (i), note that if a n = uj(F) — u n and b n = u>(F), 
then 

W^(a n x + b n ) = nF(uj(F) + x(u(F) - u n )) A 1. 

Since u n = inf{t : F(t) < n^ 1 }, we get that if h = uj{F) — u n , then as n — > oo 
we have 

F(u(F)-h)=F(u n )~n- 1 . 
It follows that if x < 0, then 

lim nF(uj{F) + x(u(F) - u n )) = ]im( F(u(F) - /i)) _1 F(lj(F) + xh) = \x\ a . 

n— *oo h—*0 'ii 

Hence 

lim F^(a n x + b n ) = \x\ a A 1 

n— >oo 

if x < 0, which proves (i). 

To prove the converse, that (i) implies (iii), we mimic the proof of the same 
statement in the Type II case, this time using the proof of Proposition 1.13 
in [14], pages 59-62. 

Assume that for some choice of a n , b n £ R, a n > 0, we have 

lim F Un (a n x + b n ) = G(x) 

n — >oo 

for i6M. Then if U = 1/F and — 1 < x < 0, we have 

lim n~ U(a n x + b n ) = \x\~ a . 

n— too 

Hence if V(y) = U*~{y) = inf{s :U(s)> y}, we infer 
\mxjy{ny)-b n )a- 1 = -y- 1 l a , 

where < y < 1 . 
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With a(t) = 0[t] and b(t) = 6[ f ], we have for < y < 1, 
lim(y(ty)-6(t))(a(t))- 1 = l-y- 1 / a 

and for < ?/2 < 2/1 < 1 , 

Km (V(t yi ) - y(t 1 / 2 ))(a(t))- 1 = y 2 " 1/a - y^. 

c — >oo 

This implies in turn that for x > 0, if < xyi < 1 and < y% < yi < 1, then 
Km a( te )(a(t))- 1 = Km ^i) ~ Hm 



t^oc i^oo a(t) t^oo F(txyi) — V{tXl)2) 

= ((xy 2 y 1/a - {x yi y l h{y2 1/a - vi l/a r l = 

and thus a(t) is regularly varying at oo of index —1/a. 

Next we show that to(F) < oo. Note that u>(F) = V(oo) = ]im y ^ 00 V(y). 
Since a(-) is regularly varying of index —1/a, we infer that if 2 _1 / Q < A < 1, 
then for some constant A > we have 

o(2 n ) < A\ n for n € N. 

On the other hand, using 

Km (V(2 n+2 yi ) - V(2 n+2 y 2 ))(a(2 n ))- 1 = y~ 1/a - y~ 1/a 

with yi = 1/2 and 2/2 = 1/4 gives that for some B > and some no € N, we 
have, for n > no, 

y(2 n+1 ) - y(2 n ) < £a(2 n+2 ) < ABX n+2 for n G N. 

Clearly, since < A < 1 , this implies 

lim V(2 n ) - V(2 n °) = Y (V(2 n+1 ) - V(2 n )) 

n>no 

< y ABX ' l+2 < °°, 

n>no 

that is, < oo. 

Since V{ty\) — V(ty2), like a(t), is regularly varying of index —1/a, we 
infer as in the Type II case that V(ty) — V(t) is also regularly varying of 
index — 1/a for all y > 0, y^l. 

With these preparations the next step is to prove that (uj(F) — Vit))" 1 
is regularly varying of index 1/a and thus that F(uj(F) — as a function 
of s is regularly varying of index a [this means F(lv(F) — u) is regularly 
varying of exponent a at 0] . The proof can now be completed using the last 
part of the proof of Proposition 1.13 in 1.3 of [14], pages 61 and 62. □ 
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7. Peaks over threshold. The probability distributions we found as pos- 
sible limits of free extremal convolutions are well known in statistics. They 
come under the name generalized Pareto distributions (see [7], Section 3.4, 
and [13], Section 1.4, for a textbook treatment of this subject). Statisticians 
have introduced a convenient parametrization of these distributions. 

Definition 7.1. For 7GR the standard generalized Pareto distribution 
is defined by its distribution function G 7 given by 

G 1 {x) = l-{l + 1 x)- lh 

for 

x > 0, if 7 > 0, 

< x < ItI -1 , if7<0, 

and if 7 = 0, Go 0*0 = 1 — e~ x for x > 0. 

These distributions appear as limits in the peaks-over-threshold approach 
to extreme-value theory, which we sketch briefly. 

Definition 7.2. If X is a random variable with distribution function F, 
for u < uj(F) the exceedance (or excess) distribution function at threshold 
u is given by 

F [u] (x)=P(X <u + x\X>u) 
F(u + x) - F(u) 



1 - F(u) 



for x > 0. 



The main result of P.O. T. theory due to Balkema and De Haan [2]: 
If F^(a u x + by) has a continuous limiting distribution function as u] 
ui(F), then 



lim 

u-[uj(F) 



FW(x)-G 1 1 



cr, 







for some shape and scale parameters 7 and a u . 

If X±, . . . , X n are n i.i.d. random variables with common distribution func- 
tion F and if > X^ > ■■■> X^ are the order statistics, that is, the 
values of the variables ordered from largest to smallest, the traditional ap- 
proach to extreme-value theory studies the distribution of the maximum 
JfW or of the first k maxima X^\ . . . ,X^ properly normalized, when 
n — > 00. The P.O.T. approach considers the distribution of the variables 
conditioned on being larger than a large threshold, which is very close to 
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the free extremal convolution studied here (see the paragraph preceding 
Definition 6.2). 

The only (technical) difference between the P.O.T. approach and the free 
probability extreme- value theory is that the latter does not fix a threshold u 
for the random variable but rather for the value of the tail of the distribution 
function. A more serious difference is that the free approach introduces a 
binary operation on distribution functions and in some sense exhibits the 
P.O.T. theory as a result of the iteration of this operation. 

The P.O.T. theory is very useful in various fields of statistics (insurance, 
reliability among others) and has been developed by R. L. Smith [15], A. C. 
Davison and R. L. Smith [6] and J. Pickands [12]. 

One should also notice that another (related) classical occurrence of the 
generalized Pareto distributions is as intensities of limiting Poisson point 
processes of extreme- value theory (see [14], page 210, Corollary 4.19). More 
precisely, F is in the free domain of attraction of G 7 iff the point measure 



converges to a Poisson point process with intensity G 7 . 
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